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Abstract 

The magnetic field generated by an electron bound in a spherically sym- 
metric potential is calculated for eigenstates of the orbital and total angular 
momentum. General expressions are presented for the current density in 
such states and the magnetic field is calculated through the vector potential, 
which is obtained from the current density by direct integration. The method 
is applied to the hydrogen atom, for which we reproduce and extend known 
results. 



I. Introduction 

Recently, Gough [1] has presented calculations of the magnetic field pro- 
duced by a hydrogen atom in various angular-momentum eigenstates. The 
calculations were done by solution of the differential equation for the various 
multipole components of the vector potential. The purpose of this paper is 
to present a method for similar calculations, where the field is calculated by 
a direct integration from the current density. The multipole expansion of the 
latter involves only the angular part of the wave functions and we shall show 
results valid for angular-momentum eigenstates irrespective of the radial part. 
Our method involves some manipulations of spherical harmonics, which may 
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appear tedious, but provide a welcome opportunity to develop skills, which 
are of use in other applications of quantum mechanics. In particular, we 
gain some insight in the systematics of occurrence of various multipole com- 
portment of the magnetic field, using identities known from the derivation of 
selection rules in spectroscopy. 

The magnetic field generated by electrons occupying states in partially field 
shells is of interest in magnetism. Magnetic neutron scattering is mainly 
due to the interaction of the neutrons's spin with the magnetic field. The 
magnetic form factor [2] is said to be given by the Fourier transform of the 
magnetization density of an atom. However, to avoid ambiguities involved 
in the definition of the magnetization density of an atom, it is important to 
bear in mind that it is the B field due to unpaired electrons that is being 
measured. Of course, the electronic states in a magnetic material are differ- 
ent from the ones bound to free atoms. Nevertheless, the study of the latter 
is useful to clarify the concepts and methods involved in the calculation of 
magnetic form factors. In view of their relevance to magnetic materials, we 
shall discuss d and / states in particular. 



II The vector potential 

The vector potential A satisfies the Poisson equation 



with the current density j as source. In reference [1], this differential equa- 
tion has been solved for the dipole and octupole components of the azimuthal 
current density. Each Cartesian component of eq. (1) is of the form familiar 
from electrostatics as one relating the potential to charge density. Accord- 
ingly, the solution can be written in the integral form 



V 2 A = - 



(i) 





where i stands for x,y and z, and R — |r — r'| . The factor can be expanded 
[3] as 




(3) 
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for r > r', where Pi is a Legendre polynomial of degree /.This equation 
is also valid for r > r after interchanging r with r . The nature of the 
angular-momentum eigenstates to be studied in the coming sections makes 
it convenient to work with spherical components. In particular, it will be 
seen that the current density has no r or 9 components and therefore, since 
divj must vanish, it does not depend on 0, i.e., j = j(r, 9)(p. Consequently, 
the Cartesian components take the form 

jx = -j(r,9) sin0; (4) 



j v = j(r,e)coa<f>. (5) 

Substituting eqs.(3) to (5) into eq.(2) then gives the Cartesian compo- 
nents of the vector potential, 



r =0 

oo 



(6) 



1=0 



Ay = ^lfff {j{r\9')cos ( t; l -Y j {^) l P l {cosa))dV' 



r =0 



00 

oo 

./ 1 , T 



1=0 

(7) 



+ // / U(r, 9') cos0% (^) l Pi(cosa))dV'} 
J J J r l=Q r 

r =r 

The function P;(coso;) can be expanded [3] as 



P,(cosa) = P,(cos0)P,(cos0') + 2 V ^— ^ PT (cos 0)PT (cos 0') ^ 

^( l + m V- (8) 

x cosm(0 — 4>) , 

where the P™ are associated Legendre functions. We shall decompose the 
current density into multipole components as 

3= £ j l (r)P^(cos9), (9) 

i=l,3,... 
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where ji(r), js{r), ... are the dipole, octupole, 32-pole,.. .component of the 
current density, which enables the use of the orthogonality of the associated 
Legendre functions, 

£p™(cose')p™(cose')sme'de' = (10) 

together with the integrals 

/ cosm(0 — ) sin d<p = 7rsin0<5 mi i (11) 
Jo 

and 

cosm(0 — ) cos0 d(f> — 7r cos <j>5 mj i . (12) 

The latter ensure that the expansion (9) is limited to m — 1. It is easily 
shown that for the dipole component the result is 

A X1 =-^sin0P 1 1 (cos^)[l f ^(rV'V + r f°° 3l {r)dr] (13) 

<J r Jr'=0 Jr'=r 

^ = ^cos^cos^M- / MrY'dr' + r / h{r)dr\ (14) 

<J r J/=0 Jr'=r 

The dipole component of the vector potential is then of the form A x = 
Ai(r)P 1 1 (cos^)0, where 

A 1 = ^-[— ji(r')r dr' +r ji(r')dr']. (15) 

<J r Jr'=0 Jr'=r 

Following the same procedure, we find the further coefficients in the multipole 
expansion of the vector potential, 

A= Mr)Pt(cos6)l (16) 

i=l,3... 

to be given by 

1 = 2T+r^H+ r y^o 3i{r )r dr+r L r — dr ] - (i7) 
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Ill Current density in eigenstates of the angu- 
lar momentum L and the spin operator S z 

III. I Multipole expansion of the orbital current density 
in angular momentum eigenstates 

The orbital current density generated by an electron in the \n,l,m > angular- 
momentum eigenstate has been given by Gough [1] as 



f = -2fi B \^ nlm \ 2 —-, (18) 
v rsmV 



where is the Bohr magneton and refers to the azimuthal component, 
the other components being zero, j° = jg = 0, and 



iJnlm = Rnl(r)Y l m (9, ( p) (19) 



is the normalized wave function . The function can be factorized into 
radial and angular parts, 



d2 rnY m 

AO _ 9,. n nl v m* mI l 

R 2 (r) 1 ' 

r 

where we have defined f lm {9) = (-l) m Ff m ^ . It is clear within the for- 
malism used in the previous section, that the various multipole components 
of the current density give rise to the corresponding components of the vector 
potential. Therefore, it will be convenient to carry out the angular integrals 
for each multipole separately To this end, in the present section we shall 
carry out the multipole expansion of the angular part of functions of the form 
(20). In doing so, we encounter the compartments 

sin0 = P 1 1 (cos0) (dipole) ; (21) 



2 

4 cos 2 9 sin 9 - sin 3 9 = -P^ (cos 9) (octupole) ; (22) 
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8 cos 4 6 sin 6 - 12 cos 2 9 sin 3 9 + sin 5 6 = — (cos 9) (32 - pole) ; (23) 

which we identify with the appropriate Legendre functions of cos 9 . To find 
the desired expansion coefficients for the angular part, first the identity [4] 



= - l^W(l-,n - m - m) e-- % T (24) 
+ ^(l + m-m + my^Y™- 1 ] 

will be used and subsequently the expansion 



y-™_v- / (2/ 1 + l)(2/ 2 + iy 
h h "f^V M2L + 1) ^ 



LO (~lLM yM /r>r\ 



where the Clebsch-Gordan coefficients obey the following rules : 

QiM 2 o = ' unless li + l 2 + L is even ; (26) 

C hm 1 i2m 2 = > wn/ess |Zi -Z 2 | < £ < Zi + Z 2 emd mi + m 2 = M. (27) 
Collecting terms, we find 

L=2«-l 



f <o\ ( ^m ( 2l + 1 ) V" r (l-m-l)(l-m) L1 

flm(9)-(-l) —^—J^ C^ 10l0 {^ j^-—^ C t _ 



1 m+1 / — m 



(Z + ro-l)(Z + m) L _, lp i, m 

!) G l-lm-ll-mri( COS ^ > 



(28) 



where we have used the defining equations 



and 



YT = (-ir[( (2 ^ ( ^ m) 7 )! ]^r(cosg)e^ (29) 



P r m (costf) = (-l) m ^p^Pr(cos9) . (30) 
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Equation (28), together with the rules (26) and (27) reveals some regular- 
ities in the occurrence of various multipole components . First, it is clear that 
only the M = 1 components appear and secondly, L — 1, 3, ..21 — 1 . This is 
in accordance with the findings of Gough, in particular that in the |2, 1, 1 > 
state only the function (21) appears, whereas in |3,2, 1 > and |3, 2,2 > the 
functions (21) and (22). It is important to note that this regularity was seen 
to follow from the angular dependence of only . This implies that the 
result will hold for any p, d, f, ect. states, irrespective of the radial function 
R(r). 

Substitution of the Clebsch-Gordan expressions leads to Table 1. 

Table 1. The coefficients of the a 1 ™ in the multipole expansion of the 
orbital current density , see equation (2.52) 









L 




I 


m 


1 


3 


5 


3 


1 


3 
8 


7 

24 


5 

24 




2 


6 

8 


7 

24 


-4 
24 




3 


9 

8 


-7 
24 


1 

24 


2 


1 


3 
8 


2 

8 






2 


6 

8 


-1 

8 




1 


1 


3 
8 







Again we find that some of the results of reference [1] are quite general: p 
states in general give rise to pure dipole fields, d states generate an octupole 
field as well, which is opposite in sign and stronger by a factor of two for 
mi = 1 than for m/ = 2, etc. These and similar regularities found in Table 2 
follow from the fact that the L-th order multipole component of the current 
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density operator is an irreducible tensor operator of rank L. 



III. II Multipole expansion of the spin current density 
for the eigenfunctions of the spin operator S z 

The current density associated with the electron spin for eigenfunctions of 
the spin operator S z has been given in reference [1] as 

fe=f r =0; 

4 = 2m^ B (sm0-^- +cos0^)W , (31) 

where 

= RljY^Y™ . (32) 
As in the previous section (III. I), eq. (25) will be used to expand Y^Y™* as 



yrnyrn* = ( _ 1)m £ (2 * + 1) C^C^Y^B, 0) . (33) 

l V / 4tt(2L + 1) 

and eq. (29) will give Y£ = (^j-^P^cos 0) where Pl = P°; consequently, 
IT^r = (-l) m ^±ll £ q&C^P^cosfl) . (34) 



To get the multipole expansion of the spin current density given in eq. 
(31), we need to express sin6Pz, and cos 0-^P L as linear combinations of 
Legendre functions. For the former, the identity [3], 

sin BPT = ^ {PR 1 " PR 1 ) (35) 
will be used, as a result, 

^0Pl = ^~ [ (P 1 l+1 -PI- 1 ). (36) 
For the latter, the definition of the M — 1 associate Legendre polynomial, 

±P L{cos e) = -pl (37) 
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will be used and subsequently the identity [5] 

cosMT = 2^pY {(' " ™ + + ( Z + ( 38 ) 

which gives 

cosOPl = ^L_{LP2 +1 + (L + 1^}, (39) 
Combining eqs. (37) and (39) we find 

m -2lTTH« + < L + 1 > F -.}- < 40 » 



cos 0- 



Substitution of eq. (34) into eq. (32) and then (32) into eq. (31) with the 
help of eqs. (36) and (40) gives 



L=0,2,... 



X 



( ® R nl LR nl \ pi ( ®Rnl , ( T , t\ &nl \ pi \ 



(41) 



By shifting the value of the summation index L, and noticing that = 0, 
the above equation can be rewritten as, 



3<t> 



r. 



L=2l+l 



v.. m i w"o (2[+l)J \^ f 1 f ./. u v ./. to ^-^n/ 

s ^^^i ^ \2L^1 10,0 

L L=l,3,.. 

2 L=2J-1 o d2 



<9r 



-(L-l)%p;- JJ ' Ctf'Cffi^QL (42) 

L=l,3,.. 

+ (L + 2)^)Pi|. 

Here again we can recognize a feature of Gough's results as being the con- 
sequence of a general rule: the highest multipole component of the magnetic 
field associated with spin currents is of order 2/ + 1. That is, s states generate 
only dipole fields, p states dipole plus octupole, ect. As we find different lin- 
ear combination of ^R n i and -f 1 in the different multipole components, no 
general expressions can be given for the ratio of different components. This 
feature is inherent in the expression (31), which shows that the spin-current 
density, unlike the orbital one, eq.(20), cannot be written as a single product 
of an r-dependent and an angle- dependent factor. 
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IV Current density in eigenstates of the total 
angular momentum 



The eigenfunctions of the total angular momentum being superpositions of 
the two eigenfunctions of S z , the calculation of the expectation values of the 
current densities involves two-component spinors. As in the previous case, 
we find that the spin current density is not factorized in radial and angular 
parts. However, it turns out that the total current density can be written as 
a single product of r- and ^-dependent functions. Therefore, we shall derive 
this expression first, in section IV. I, and turn to the problem of the multipole 
expansion of the radial part in section IV. II 



IV. I The factorization of the total current density 

The eigenfunctions of the total angular momentum are of the form 



where 



^(j = l-±, mj = m+±) = Rni^tYr, 
M = /-|,m,=m+i) = -R nl sJ^Lrr +l 



for the low-lying spin-orbit coupled state and 



M = l + \,m 3 =m + \) = R nl ^^Yr, 
M = l-i, mj = m+i) = Rni\f^Yr +1 



(44) 



(45) 



for higher-lying state. The current density operators applicable to wavefunc- 
tions of this form are represented as 2 x 2 matrices: 



)° = — - ttVtt*) (46) 



where V stands for I ^ ^ ) , and 



V 

j s = -^ B Vx < a >, (47) 
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where a stands for the Pauli matrices, 



I l)> a v = ( i o) ; ^=( J -J)' (48) 

which can be transformed to spherical components using the transformation 
matrix 

sin 9 cos sin 9 sin <fi cos 9 

cos 9 cos <fi cos 6> sin — sin# | | cr y ] (49) 
— sin cos 6 



giving 



0> 



ay > 






1- 






/ cos# 


y sin#e' 


-( 


' 







sin 9e %<i> \ ( — sin 9 cos #e 



-cos# J ,ae \ cosfle 1 * sin# 



-i e ~i<t> 





(50) 



As a consequence of the diagonal form of the operator (46), the orbital current 
density is found to be the sum of the currents generated by the separate spin 
components, given by eq. (20): 



j° = -^W T >T + (m + 1)^)0. (51) 



(52) 



Substitution of the spinor components (44) gives 

+ (/ + m + l)|^| 2 ^}0. 

The evaluation of the spin current density is a bit more elaborate as it requires 
the three components of the spin density . Taking the expectation values of 
the matrices (50) using again the spinor components (44), we find 



<a r > = R 2 nl ^p [ { cos0[(Z - m)\Yr\ 2 — (l + m + 1)|YT + T] 
-2sm6y/(l-m)(l + m + l^Y^+V^}; 

<a e > = R 2 nl ^{ ~ an6[(l - m)\Yr\ 2 - (I + m + 1)|YT + T] 
- 2 cos 6y/(l -m){l + m + l)^ m *F z m+1 e -^}; 



(53) 



(54) 
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< (70 >= 0, (55) 

where we have made use of the fact that y™*y™+V^ = Y l m Yj n+1 *e i4 ' is 
real. As < > vanishes and < a r > and < a e > are independent of , 
eq.(47) will provide only a component, which depends only on r and 9 : 

r = - < °e » ~ | < * r >}l (56) 



To evaluate this quantity we require the derivative 

d_ 



-(l + m+ l)^|^ m+ T] - 2 sin 6 y/ (I -m)(l + m + l) (57) 
Substituting this expression and eq.(54) into eq.(56), we find 

r , 1 = _ ^ f r ggli(_ sin gr ( /_ m) |^|2 

J ^=^2 r (2/ + 1)1 * I LV ;| ' 1 

- (l + m + l)|^ m+1 | 2 ] - 2 cos 9^{l - m )(l + m+ VjY^Y^ 1 e'^ 

- Rl\ cos9[{l - m)^\Yr\ 2 -{l + m+ 1)^|>T + T] 
-2sm6y/(l -m)(l + m + 1)— ^ m *^ m+1 e~^})0. 

(58) 

To eliminate the derivatives from the second pair of braces, we make repeated 
use of the identities [4] 

Ay™ = m C ot 6Y™ + y/(l -m)(l + m + l)^ m+1 e -^ 

cw (59) 

= -m cot 0*7" - v /(/ + m)(/-m + l)r / m - 1 e ie 

and we find that the awkward terms containing in eq.(52) can be com- 
bined with terms, so that the seemingly disparate angular functions 
appearing in j° and j s provide the desired factorized from for j = j° + j s : 

- (/ + m + l)|^ m+1 | 2 ] + 2 cos 9^ (I - m)(l + m + l)F™*F™ +1 e -^}0. 

(60) 
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This indeed is a product of a radial and an angular function, the former being 
independent of m. Similarly, for the case j — I + \ we find 



h=«\ = wfi [ ^ ' 2l ^ ] i shl9[{l +m+ mr? - (/ " m) ^ m+1 ^ 



5 21 + 

- 2 cosO^/il -m)(l + m+ l)^ m *^ m+1 e -^}0. 
IV. II Multipole expansion of angular part. 



(61) 



To find the multipole components of the total current density, as given in 
eqs.(60) and (61), we again use eqs.(25) and (29). These readily yield eq.(34) 
and 



:-i) 



m+1 



j~ J n L r<L 1 pi 
PFTr — I — TT '0 i0°! -m / m+l r L- 



4vr ^ v/L(L + l) 
To come to our final result, we require eqs.(36) and (39) and the identity [4] 

y/L(L + l)Ci} m l m+1 = a/(/ — m){l + m + 1) (Cfl°( m+1 ) ; m+1 + Cf_° m j m ) . 

(62) 

As before, some shifting of summation indices is necessary, in order to group 
terms according to the order of the Legendre function they belong to and 
this results in: 



1) 



47T 



<9r 



L=l,3,.. v 7 



2(Z + m + l)\Cf-\ °_ m + (Z + m + 1)[L - 2(Z - m)]C</^ _ (m+1) }^ 



2/-1 



E 

L=l,3,. 



(L+l)(2L + 3) 



CftiW - m)[L + 1 + 2(Z + m + l)]C7/£ °_ T 



+ (Z + m + 1)[L + 1 + 2(Z - m)]C^l° u _ (m+1) }Pi 



(63) 
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With similar manipulations, we get from eq.(61) 

21+1 



[E 



1 



(-l) m \ dR 2 nl Rli 
^ 4tt L dr r J L ^ L(2L - 1) 

L— 1,3,.. 



+ 2(Z - m)]C^ m ) ° m + (Z - m)[L + 2(1 + m + 1)]C^ _ (m+1) }Pi 



2J-1 

E 

L=l,3,. 



1 



(L + l)(2L + 3) 



Cffi °{(l + m + 1)[L + 1 - 2(1 - m)]C%) °_„ 



+ (Z - m)[L + 1 - 2(Z + m + 1)]C/^ _ (m+1) }Pi 



Equations (63) and (64) can be further reduced to : 

2/-1 



(64) 



ij=l-^,mj 



■1) 



47T 



<9r 



+ 2 « + 1 >T 1 ] E [l(^I) C -o°{('— )[L 



L=l,3, 



2(Z + m + l)]C7f-i ° m + (Z + m + 1)[L - 2(1 - m)]C/^ _ (m+1) } 
r^o + /o°{ (* - m) [L + 1 + 2(Z + m + 1)]C?£ m 



(L + l)(2L + 3) 
+ (Z + m + 1)[L + 1 + 2(1 - m)]C^ }1 P'J 



(65) 



hi) 

An 



21-1 



m r Qp2 p2 
UJX nl _ 21 — BL 

dr r \V L ^\L(2L-l) 



E 



+ 2(Z - m)]Cfa) °_ m + (l-m)[L + 2(l + m + 1)]C/^ _ (m+1) } 
CftzVH (' + m + 1) [L + 1 - 2(Z - m)]<7/£ ^ 



+ 



(L+ l)(2L + 3) 
1 

(2/ + l) 2 



+ (Z - m)[L + 1 - 2(Z + m + 1)]C, L+1 



(4Z + l)(Z + m + l) 



i m+l Z -(m+1) 

r 2l ^2/ pi ll 
°Z J 0°Z m I -m r 2l+l <Pi 



})Pl 



(66) 



where we have make use of the identity [4] (Z + m + l)Cf 



z m+l Z —(m+l) 
a+6 a+/3 



0bf3 



to 



(Z — mjCi l^i _ m , which follows from the general expression for C a a 
eliminate and reduce the highest-order term sum in eqs.(63) and (64). These 
two equations may look discouraging but eminently suitable for numerical 
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calculation of particular cases, as will be demonstrated in section V. Also, 
the coefficients of the multipole expansion are universal in the sense the ones 
given in table 1 were found to be. 



Table 2 The coefficients a 3 ™ 3 in the multipole expansion of the total cur- 
rent density, see equation (2.54). 
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35 


15 
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Table 2 lists such coefficients according to the definition 

Wi, = ±£ + P t W + Dl^} E <M «*> 

L=l,3... 

The notation in this expression implies what is not immediately obvious from 
eqs. (65) and (66): the expansion coefficients are independent of I. In other 
words: for instances, no matter whether a j = 5/2 sextet involves a d or 
/ electron, each of the three allowed multipole components of the current 
density will have the same relative strength for different rrij values in both 
cases. Again, this is a consequence of the Wigner-Eckart theorem, which 
is applicable, because the total current density has the same property as 
the orbital one: its multipole component of order L is an irreducible tensor 
operator of rank L. 

V Examples 

V.I The magnetic field generated by the orbital current 
in the state \n = 3, I = 2, mi = 1 > 

Since the radial part of the wave function of the state |3, 2, 1 > is 

R 3 2(r) = ^ e^r 2 , (68) 
1215a/ % 

after substitution of R 2 2 and /2i(#) from table 1 into eq.(20) we get, 

H 5(3 9 )(a£) V 8tt 1 8tt 3J K ' 

Comparing the above equation with eq.(9) gives, 

•?i( r ) = 5( 3*xV e ^ r3 and j3(r) = WM> e ^ 

into eq.(17) and integrating using the identity [1] 

J r n e~ r/a dr = -ae~ r/a {r n + nar n ~ x + n(n - l)a 2 r n - 2 + ... + n\a n } (70) 
we have , 



/:l ( r) = 5(3*)(4)n e3a ° r3 and M r ) = 5(3 9)( 4 a v )7r e 3a or 3 . Substitution these js 



_ fipfiB sin 9 ( =*r, 2 r 3 S 

ATra 2 lea ° l (5)(36)a 3 + (5)(35) a 2 

19 r 2 2a a 2 ,. a 2 ,, 

J 1 1 ^ _L _U) U 1 

(5)(3 4 )a 9 3 r r 2j r 2 J 



(71) 
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27/!o/Ub(4 cos 2 9 sin 9 — sin 3 9) < =2l 1 r 3 4 r 2 
3 = 2^ |6 "° ( (5pj ^ + W&)$ m 

4 r 2 4 g 2 a 2 2 a 3 a 4 a 4 
(5)(3 6 )a 3 5 3 4 r 9 r 2 3 r 3 r 4 ^ r 4 J 

The magnetic filed B = curl A has two components [1], B r = rsi 1 n6 , ^(sin 9A^), 
and ^ = ~-^{rA^). Evaluation of these components for the ^ and A 3 re- 
sults in the magnetic fields , B ri , B ra , B dl and Bg 3 where, 



_ wb cos 9 ^ 2 r 2 



2vra 3 V (5)(3 6 )a 2 (5)(3 5 )a 

!9 , 2 ao 2a| a| a|, 
(5)(3 4 ) 9 r 3r 2 r sJ r 3 ] 



(73) 



5 



54/i 0/ ub(5 cos 3 9 



rz 



o 



10) , ^ 
— - e 3a o 



(5) (39) a 2 (5) (3») do 



+ 



(5)(3«) 



_2_ao _4 a| 

35 r 34 r 2 



2 of 2a| 
9 r 3 3 r 4 



_2 \ 

.5 J 



(74) 



/i /x B sin 6 1 f 4 r 3 8 r 2 



47ra 3 1 (5)(3 7 )a 3 (5)(3 e 



2 



14 r 22 2a 2 a 2 a 3 a 3 
j 1 1 _| 1 _|_ " I 

(5)(3 5 )a 5(3 4 ) 9 r 3 r 2 r 3 ^ r 3 J 



81u /is(4cos 2 6'sin6' - sin 3 6») , ^r. 2 r 3 4 r 2 
= ^—3 i e3ao (t^TTTTs + 



27rag L v (5)(3 n )ag (5) (3 



10 ^a 2 



(76) 

2 r , 4 , 2_ao £^0 , 2 a| 2 a| a| a|, 
(5)(3 8 )a (5)(3 6 ) 3 5 r 3 4 r 2 9 r 3 3 r 4 r sJ r 5 1 



This example serves as a verification that our method reproduces the results 
of ref.l, which it does. Figure 1 illustrates the result in terms of field lines. 

IV. II The magnetic field generated by the state \n = 

3,/ = 2, j = 3/2, rrij = 3/2 > 

Substitution of the radial wave function R 32 of the previous example and the 
coefficients a 3 ™ 3 of table 2 into eq.(67) gives 



17 



-2r 

J M = ^m4 { 15a °- r)( - Pl+ 35 Ps) 

Following the same procedure of the above example, the expressions for 
vector potential and for the magnetic are, 

3fi fi B sm9 f 4 r 4 2 r 2 o - o 



lOvra 2 (5) (3 s ) a 4 (5)(3 4 )a 2 (5)(3 2 ) a 9 



-27/i /is(4cos 2 6'sin6' - sin 3 6>) , =*l 4 r 



2 



20vra 2 (5)(3 8 )ag (5)(3 6 )a 

+ — + 1^ + H^l + H^l + ^ _ 2 r \ _ M\ 

3 5 3 4 r 9r 2 3 r 3 r 4 (5)(3 10 )ao r 4 ' 

g^BCQsg r ^ ( 2 r 2 2a 2a 2 

IOtto 3 1 1 (5)(3 4 ) a (5)(3 2 ) 9 r 3 r 2 

r 3 (5)(3 8 )< r 3 J 

-27/^0^(5 cos 3 6» - 3cos0) r =2r, 4 r 4 
^3 = l e3a ° 



5vra 3 V (5)(3 8 )a (5)(3 6 ) 

+ 4. 1^1 4. H^l 4. H^l 4. 4 _ 2 r \ _ 4\ 

3 5 r 3 4 r 2 9 r 3 3 r 4 r 5 (5)(3 1( Vo r 5 1 

3/io/i_B sin j- ^2r 4 r 3 4 r 2 2 r 

+ _A_ 4 2 ^ 4 2 ^I 4 4 _ 8 A _ ^0 1 

5(3 3 ) 9r 3r 2 r 3 (5)(3 9 )a 4J r 3 J 



81/x yM4cos 2 6>sin6> - sin 3 6>) , ^ 2 r 



-D(9i = 5 i e 3a ( — — t + 



207ra 3 L v 3 n a|] (5)(3 10 )a 2 

4 r 4 2 a 4 a 2 2 a 3 , 2 a 4 
+ (5X38) ^(5)(36) + 3^7 + 77 + 9^ + 37 + 7 
4 r 4 a| 

> r 5 i 



(5)(3 12 )a 
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Figure 2 is a graphical representation of the magnetic field in this case. 
The pattern is much simpler than the previous one; it is reminiscent of the 
field of a simple current loop. The same holds for the case j — § rrij — §, 
as seen in fig. 3 . This becomes understandable if one notices, using the 
coefficients of table 2, that j(r,9) is proportional to sin 3 # for j = rrij = | 
and to sin 5 9 for j = rrij = |. Thus, with increasing j, the current tends to be 
confined in the 'equatorial' plane. Comparing figs. 3 and 4, one realizes the 
limits of the 'universality' of the relations summarized in table 2. Although in 
these two cases the ^-dependence of j is identical, the different r-dependences 
influence the resulting pattern of field lines in a qualitative way. 
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